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Abstract 

We present a new rigorous method for estimating statistical quan¬ 
tities in fluid dynamics such as the (average) energy dissipation rate 
directly from the equations of motion. The method is tested on shear 
flow, channel flow, Rayleigh-Benard convection and porous medium 
convection. 


1 Introduction 

One of the most fascinating features in turbulent flows is the emergence of 
complicated chaotic structures involving a wide range of length scales be¬ 
hind which “typical” flow patterns are still recognizable. The state of motion 
is too complex to allow for a detailed description of the fluid velocity and 
experimental or numerical measurements of certain system quantities ap¬ 
pear disorganized and unpredictable. Yet, some statistical properties are 
reproducible [9]. One of the challenges in theoretical fluid dynamics is thus 
the derivation of quantitative statements on turbulent flows. Many of the 
approaches consist of various approximation procedures, the imposition of 
physically motivated but ad hoc assumptions (e.g., “closure”) or the intro¬ 
duction of scaling hypotheses. Rigorous results beginning with the equations 
of motion directly are therefore indispensable for checking the validity of the 
imposed simplifications and for justifying secondary models and theories. 

The idea of extracting information about driven turbulent flows via bounds 
on physical quantities through mathematically justifiable operations and 
without imposing ad hoc assumptions goes back to the pioneering works 
of Malkus, Howard and Busse in the 1950s and 60s PI El an Slug. These 
authors applied variational approaches for the derivation of bounds on the 
energy dissipation rate in models for shear flow and heat convection. In the 
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1990s, Constantin and Doering introduced a practical framework for estimat¬ 
ing physical quantities rigorously and directly from the equations of motion, 
which they called the “background flow method” 

The background flow method is an extremely robust method for con¬ 
structing bounds in fluid dynamics and builds up on techniques developed 
by Hopf to generalize Leray solutions of the Navier-Stokes equations to fi¬ 
nite geometries with physical boundary conditions [12]. In this method one 
manipulates the equations of motion relative to a steady trial background 
state. Decomposing the quantity of interest, e.g., the energy dissipation, 
into background profile, which satisfies the forcing conditions, and fluctuat¬ 
ing component, the background part yields an upper bound if the fluctuation 
term satisfies a certain nonnegativity condition, which is often referred to as 
the spectral constraint. In a certain sense, finding the least upper bound 
using Constantin and Doering’s method resembles a variational saddle point 
problem. The equations imposed on the fluctuations necessarily include the 
equations of motion. To simplify the derivation of the Euler-Lagrange equa¬ 
tions corresponding to the variational problem for the fluctuation term, how¬ 
ever, fluctuations are often chosen in a much larger class of functions. In 
other words, the spectral constraint in the background method is required 
to hold for an infinite-dimensional set of vector fields, that strictly contains 
the solutions of the equations of motion. In that case, enforcing the spectral 
constraint may possibly yield to an overestimation of the quantity of primal 
interest. 

After its introduction, the background flow method started its triumph 
as the upper bound method with applications ranging from various prob¬ 
lems in turbulent heat convection and boundary-force and body-force driven 
turbulence to idealized models in magnetohydrodynamics. The method was 
soon improved by Nicodemus et al. who introduced an additional balance 
parameter [19], and Kerswcll showed that this improved method is actually 
“equivalent” to the approaches of Busse and Howard [16] . 

Apart from its practical performance, for many years, the background flow 
method was considered as a rigorous manifestation of Malkus’s marginally 
stable boundary layer theory [18]. The latter is based on the assumption that 
turbulent boundary driven flows organize themselves into marginally stable 
configurations. If the well-mixed core is bounded by thin laminar boundary 
layers, the thickness of these layers is determined by the condition of marginal 
stability. The association of the background flow method with Malkus’s the¬ 
ory relies on the surprising (?) observation that the spectral condition in 
the backgroTind flow method resembles a nonlinear stability condition on the 
background flow. A recent work of Nobili and Otto [20], however, proves the 
failure of this association — at least in the context of infinite Prandtl number 
Rayleigh-Benard convection: The authors compute the least upper bound on 
the Nusselt number (the quantity of interest in Rayleigh-Benard convection) 
within the framework of the background flow method. This bound however, 
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exceeds the bound derived by Otto and the author in |2T] using completely 
different methods. In the context of Raylcigh-Benard convection, it thus 
seems that a physical interpretation of the background flow method is mis¬ 
leading. Whether the background flow method indeed gives physically rele¬ 
vant information (apart from scaling bounds) in different problems of fluid 
dynamics can only be speculated. 

In this paper, we focus on the energy dissipation rate as an example of one 
specific physical quantity and present a new method for its rigorous 3 estima¬ 
tion directly from the equations of motion. (In fact, the method was already 
introduced in eh by Otto and the author, but its universality was not seen 
at that time.) To allow for a straight comparison with the background flow 
method, the method is tested on the problems considered by Constantin and 
Doering in [HI j5[ 17, [8j. More precisely, we study the classical fluid dynamics 
problems shear flow , channel flow, Rayleigh-Benard convection and porous 
medium convection. All of these problems have in common that they can 
be considered as model problems for boundary-force or body-force driven 
flows or for thermal convection. In fact, we will recover the same results as 
Constantin and Doering in the above mentioned papers. 

Our new approach presented in this paper is entirely different from the 
background flow method in that it is rather based on (local) conservation 
laws. More precisely, at the heart of our method conservation laws for the 
certain flux components are established. The conservation laws are local in 
the direction of the symmetry axis but global in all other spatial and time 
variables. Typically, the conserved quantity can be explicitly related to the 
energy dissipation rate. The former can be averaged over boundary layers 
along the rigid domain walls where the velocity held is small due to no¬ 
slip boundary conditions and it can thus be easily controlled by the viscous 
dissipation rate. Optimizing the width of the boundary layer then yields a 
bound on the energy dissipation rate. How the new method applies to fluids 
with different boundary conditions, e.g., no-stress, has to be seen. 

The author believes that the strength of the new method relies on the fact 
that it uses the equations motions and some secondary derived physical laws 
directly instead of working with a rigid “upper bound construction” which 
in some cases is too restrictive to yield the optimal result, cf. [211 ED]. The 
involved mathematical operations are elementary. 

We apologize in advance that in this paper we will not apply any effort in 
explicitly computing nor optimizing the prefactors in our bounds. The aim 
of this paper is solely to advertise an alternative method to the most widely 
used background flow method to derive upper bounds on physical quantities 

*We caution the reader that our results are still only “formally” true in the sense 
that many of the performed mathematical operations apply only to sufficiently smooth 
solutions of the Navier-Stokes equations. Of course, the results can be made rigorous if the 
analysis is performed on suitable week solutions, e.g., Leray solutions, in the appropriate 
framework. 
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in turbulent fluid flows. 

The article is organized as follows: After fixing the notation in Section [2j 
we will derive upper bounds on the energy dissipation for shear flows (Section 
[3]), channel flows (Section 0]) , Rayleigh-Benard convection (Section EJ) and 
porous medium convection (Section [6]). 


2 Notation 


In what follows, we will be try to be as consistent as possible with regard to 
the notation even though different physical problems will be considered. 

Our models are nondimensionalized. 

Our system is a layer of fluid in the box [0, L) d ~ l x [0,1] where L is an 
arbitrary positive number that will not enter in our analysis. Throughout the 
article, we will refer to the first d— 1 coordinates as the horizontal and the last 
coordinate as the vertical one. We write x = (y, z ) £ R 6 ^ 1 x R accordingly, 
and denote by {ei,..., e^-i, &d\ the canonical basis of R (/ = R d_1 x R. 

We assume periodic boundary conditions in the horizontal directions for 
all variables involved. The horizontal boundaries are rigid and the imposed 
conditions will depend on the particular physical problem under considera¬ 
tion. 

The fluid velocity is denoted by u, and we write u = (v, w) £ R d_1 x R to 
distinguish the horizontal velocity vector from the vertical component. The 
hydrodynamic pressure is p and T is the temperature field. 

We consider the rate of energy dissipation 

e = [ (| Vw| 2 ) dz, 

Jo 

where (•) denotes the horizontal and time average, that is, 


(/) 


lim — [ 

rfoo T J 0 L d 1 



f{t,y) dydt. 


In general, the long-time average need not to exists, even if finite time av¬ 
erages are bounded, and we could be more careful at this point by choosing 
the lim sup instead of the lim. However, for the sake of a clearer statement 
and to simplify the subsequent analysis, we will be quite formal in most of 
our computations. 


3 Shear flow 

As the simplest example for a boundary-driven flow, we consider a fluid 
which is confined between two parallel plates that are moving at a constant 
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speed relative to each other. The equations of motion are the Navier-Stokes 
equations in a box [0, L) d ~ l x [0,1], 


d t u + u ■ Vu — A u + Vp = 0, 
V ■ u = 0, 


( 1 ) 

( 2 ) 


and, as stated in the previous section, the velocity u and the pressure p are 
both periodic in the horizontal variables. At the horizontal plates, we assume 
no-slip boundary conditions for the velocity field. If the upper boundary plate 
is moving with constant speed Re in direction e\ while the lower plate is at 
rest, we must have 


u = 0 for z = 0 and u = (Re)e i for z — 1. 


(3) 


Notice that all quantities are nondimensionalized and Re is the Reynolds 
number. 

The scaling of the energy dissipation rate e as a function of the Reynolds 
number is of fundamental importance in many engineering applications, since 
for steady states, the energy dissipation rate measures the rate at which work 
must be done by an agent to keep the upper plate moving. 

The interest in understanding the energy dissipation rate in boundary 
driven fluid flows as a function of the Reynolds number dates back to Stokes. 
It is well-known PH US] that solutions to the Stokes equation minimize the 
dissipation rate among all divergence-free vector fields with a fixed velocity at 
the boundary. In many situations of physical interest, this solution is laminar 
and continues thus to exist as a solution of the Navier-Stokes equation. In 
our case, the laminar solution is the so-called Couette flow uc{z) = (Re)ze i 
which dissipates energy £c = (Re) 2 . For small Reynolds numbers, uc is 
a stable solution of the Navier-Stokes equation in the sense that any suffi¬ 
ciently small perturbation will decrease in time. For large Reynolds numbers, 
however, perturbations of the laminar steady state are unstable and the flows 
can become chaotic or turbulent. The energy dissipation rate is a monotone 
function of the Reynolds number and acts as a measure of how turbulent 
the flow is. With regard to the presence of steady regular solutions also in 
the high Reynolds number regime—let us call these solutions “ungeneric”— 
absolute lower bounds on the energy dissipation will be dictated by untur- 
bulent flows: £ > £c = (Re) 2 . A rigorous scaling theory for the dissipation 
rate can hence only be an upper bound theory. An upper bound, however, 
sets limits on the possible turbulent structures of the flow and is thus an 
indispensable knowledge in the study of turbulence. 

Developing a conventional statistical turbulence theory for high Reynolds 
numbers, Constantin and Doering predict the “logarithmic friction law” 


(Re) 3 


as Re 1, 


(log Re) 2 
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cf. |g Appendix A], which is in accordance with the experimental data de¬ 
rived in fra. In the same paper, the authors derive a first upper bound 
which proves the conjectured rate up to the logarithmic factor, that is, they 
prove that £ < (Re) 3 using the background flow method. In the following, 
we reproduce Constantin and Doering’s bound using our new method. 

As a starting point of our analysis, we recall that the energy dissipation 
can equally be expressed, for instance, as the trace of the vertical derivative 
of Vi — u ■ ei at the top plate, that is, 


£ = Re{d z \ z=1 V!). (4) 

The identity follows from testing (JTj) with u, integrating by parts and using 
(J2]) and (J3]). Moreover, multiplying (JT]) by e±, taking the horizontal and time 
average, and using (J2]) and the periodicity, we see that 

( wv\ — d z V\) = const for all z G [0,1], 

and thus, in view of ([3]) and (J4J) , we obtain the formula 

e = Re(d z v i — wvi) for all z G [0,1]. 


In particular, averaging this identity over some boundary layer [0,£], where 
0 < l < 1 has to be determined later, we see that 


£ ( v i\ z =e) ~~j\ (' wv i > dz ’ 


where we have used the fundamental theorem together with the no-slip 
boundary condition (J3]). We now apply the Cauchy-Schwarz and Poincare 
inequalities and deduce 


£ < ^ sup (H : -) 1/2 + ^ f ( \u\ 2 )dz 

zG[0,f] t Jo 

u / ri \ 1 / 2 rZ 

< W2 ( \d z u\ 2 ) dz) + (Re)£ ^ ( \d z u \ 2 ) dz 

< ^£ 1/2 + ( R e)te. 

Optimizing the last expression in £ yields i ~ e" 1 / 3 <C 1, and thus e < 
( Re)e 2 / 3 , which entails 

£ < (^) 3 - 

This estimate agrees with the bound derived in [6]. 
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4 Channel flow 


As the simplest example for a body-force driven flow, we consider a fluid in 
a rectangular domain, which is driven by a pressure gradient in a direction 
of one of the horizontal boundary plates. The problem is modelled by the 
forced Navier-Stokes equation 

d t u + u ■ V« — Au + Vp = (G'r)ei, (5) 

V • u — 0, (6) 

in [0 ^L) d l x [0,1], supplemented with no-slip boundary conditions 

u = 0 for z E {0,1}. (7) 

Here, Gr denotes the nondimensional Grashof number. We remark that when 
the force is specified a priori as in the present situation, the Reynold number 
is an emergent quantity. 

The channel flow problem has the laminar solution Up{z) = |( Gr){z 2 — 
z)e i, the so-called Poiseuille flow, for which the energy dissipation rate is £p = 
{Gr) 2 / 12. The Poiseuille flow plays a similar role in the channel flow problem 
as the Couette flow in the shear flow problem: up solves the corresponding 
Stokes equation and is an unstable solution of the Navier-Stokes equation in 
the high Grashof number regime. When bounds are expressed in terms of the 
Grashof number instead of the Reynolds number, the Stokes limit represents 
an upper bound on the energy dissipation rate: £ < £p = {Gr) 2 / 12. 

The scaling of the energy dissipation rate in the high Grashof number 
regime is expected to obey the modified logarithmic friction law 

(Gr) 3 / 2 
(log Gr) 2 

- similar to the shear flow. Due to the presence of ungeneric solutions of 
the Navier-Stokes equation, we can only expect to derive a one-sided version 
of this scaling law. A first rigorous lower bound on the energy dissipation 
rate was established by Constantin and Doering [5], e > (Gr) 3 / 2 , which is 
optimal up to the logarithm. In the following, we will recover this bound 
following our new method. For further improvements we refer to ra and 
references therein. 

We see from testing (J5]) with u, integrating by parts and using (J6]) and (JTJ) 
that the energy dissipation rate is proportional to the average flow velocity 

£ = Gr f ( v\) dz. (8) 

Jo 

Moreover, multiplying (J5|) by e±, taking the horizontal and time average, and 
using ((6j) and periodicity, we obtain 

d z (wv i — d z vi) = Gr for all z e [0,1]. 
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From the no-slip boundary conditions (J7J) we thus infer that 


(wv i - d z v\) = ( Gr)z - ( d z \ z=Q vi) for all z e [0,1]. (9) 

On the one hand, averaging in z over the strip [1 — £, 1] of width yields 


Gr 


(v\w d z vi) dz T (d z \ z=0 o)- 


'!-£ 


On the other hand, averaging fl9]) over [0,£] yields 


(d z \ z=0 vi) ~ - / (d z vi-viw)dz + £Gr. 


Since t 1, a combination of the previous two estimates gives 


Gr<- (d z v 1 -v 1 w)dz+- {v x w - d z Vi) dz. 

"Jo " Ji-e 

We can apply the same arguments as in the shear flow case considered in the 
previous section to deduce 


Gr < 


£ l/2 


+ le. 


The optimal i is of size i ~ e 1//3 , which implies that 

(Grf 2 <e. 

In view of (l8j) . this bound is equivalent to 

(Gr) 1/2 < f\vi) dz, 

Jo 

which agrees with the bound derived in [5] and appears to be sharp to within 
logarithms. 


5 Rayleigh-Benard convection 


Raylcigh-Benard convection is the transport of heat by thermal convection in 
a fluid layer that is heated from below and cooled from above. The problem 
is modelled by the equations of the Boussinesq approximation 


d t T + u • VT - AT 

= o, 

(10) 

V • u 

= 0, 

(11) 

(d t u + u ■ X7u) — A u + Vp 

= RaTed ■ 

(12) 



The system is nondimensionalized and admits two controlling parameters, 
the Rayleigh number Ra and the Prandtl number Pr. The equations are 
complemented by the boundary conditions 

T = 1 on z — 0, and T = 0 on z — 1, (13) 

representing heating at the bottom and cooling at the top, and no-slip bound¬ 
ary conditions for the velocity field, 

u = 0 on z <G {0,1}. (14) 

The quantity of interest in this model is the so-called Nusselt number, a 
measure for the average upward heat flux. It is defined by 

Nu = [ ((uT - VT) • e d ) dz = [ ( wT - d z T) dz. 

Jo Jo 

The scaling of the Nusselt number in terms of Ra and Pr is a problem 
of enormous experimental, numerical and theoretical research for over fifty 
years. For a recent review, we refer to [T] and references therein. For all 
values of Ra and Pr , a laminar solution is given by u c = 0 and T c = 1 — z, 
which corresponds to pure conduction. The corresponding Nusselt number 
is Nu c = 1, and the laminar solution is unstable for large Rayleigh numbers. 

In the high Rayleigh number regime, Ra 1, the scaling of the Nusselt 
number is proportional to the scaling of the energy dissipation rate. Indeed, 
testing (fT2l) with u, using the imcompressibility assumption (11111 and invoking 
the boundary conditions (1T3l) and ()T4li for T and u yields that 

e = Ra(Nu — 1) ~ RaNu. (15) 

The bound on the Nusselt number in the ultimate turbulent regime is ex¬ 
pected to be 

Nu ~ Ra 1/2 , 

if Ra 1, uniformly in Pr\_ j] In the following, we derive an upper bound on 
this scaling with the help of different representations of the Nusselt number, 
similar to the approach in the previous two sections. Averaging the heat 
equation (jTOl) and using periodicity and dill) , it follows that the heat flux is 
constant on every horizontal slice, that is 

Nu = (Tw — d z T) for all z <G [0,1]. 

Rn fact, this scaling is in conflict with the rigorous bound derived by Whitehead and 
Doering for 2D Rayleigh-Benard convection with free-stress boundary conditions, in which 
case Nu < Ra 5 ' 12 could be proved [25] ■ The proof makes heavily use of the 2D struc¬ 
ture (no vortex stretching) and the free-slip boundary conditions (homogeneous vorticity 
boundary conditions). Whether such a bound can be extended to our problem at hand 
is not clear to the author. Recent numerical simulations at least indicate that the 5/12 
scaling should be expected for any solution of m d with finite energy dissipation rate 

m- 
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Now, averaging this identity over a boundary layer of thickness £ G [0,1], 
and using the maximum principle on the temperature max |T| < 1, which is 
enforced by the boundary conditions (if not initially, then exponentially fast 
in time), we obtain that 

If 1 11 f l 1 

Nu<- JI (wT) dz + - < - II ( \w\)dz +-. 

We use Poincare’s and Holder’s inequalities and (H4l) to bound the integral 
over the vertical velocity component by the energy dissipation rate e, that is 

/ i i*i / \ 

(\w\) dz < £ J (\d z w\) dz < £ 3 / 2 (J ((d z w) 2 ) dz) < £ 3 ^ 2 (RaNu) 1 ^ 2 
by (lT5li . so that 

Nu < £ 1/2 (RaNu) 1/2 + j. 

Optimizing in i yields that £ ~ (RaNu)^ 1 ^ 3 , which entails that 


Nu < Ra 1 ' 2 . 


This is precisely the same scaling law derived by Constantin and Doering in 

13 . 

Applying the same method combined with sophisticated maximal regu¬ 
larity arguments, Choffrut, Nobili and Otto [4] recently obtained new bounds 
on Nu which improve this bound in certain Ra-Pr regimes. The results in 
particular apply to the large Prandtl number regime. The Nusselt number 
bound can be interpreted as a bound on the average temperature gradient, cf. 
(l2TTl below. Developing techniques similar to those presented in this paper, 
the author derived bounds on higher order derivatives of the temperature 
held in infinite Prandtl number convection and estimated deviations of the 
average vertical temperature profile from linearity [23]. 


6 Porous medium convection 

We finally consider thermal convection in a porous medium. In this case, 
Darcy’s law approximates the Navier-Stokes equations, and the Raylcigh- 


Benard system dTOli (IT2|) reduces to 

d t T + u-\/T - AT = 0, (16) 

V-u = 0, (17) 

u + Np = ( Ra)Ted■ (18) 

As before, the non dimensional number Ra is the Rayleigh number. The 

boundary condition satisfied by the fluid velocity are 

w = 0 on z e {0,1}, (19) 
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and we suppose that the container is cooled from above and heated from 
below, modeled by 

T = 1 on z — 0, and T — 1 on z — 1. (20) 

Again, we assume periodic boundary conditions in all horizontal directions 
for all quantities involved. 

As in the case of classical Rayleigh-Benard convection, the quantity of 
interest here is the Nusselt number 

Nu — f (wT — d z T) dz. 

Jo 

Before estimating the energy dissipation rate in this example, we start with 
the study of the Nusselt number. Experiments and numerics suggest that 

Nu Ra as Ra 1, 

cf. mi. Because of the existence of ungeneric laminar solutions, we can only 
expect to prove the upper bound Nu < Ra. We first need to establish some 
alternative identifies for Nu. We Erst recall that averaging (fTfilh we obtain 
(wT — d z T) = const , and thus 

Nu = (wT — d Z T) for all z e [0,1]. 

In particular, Nu = —(d z \ 0 T). Testing now (TT6|) with T, integrating by 

parts yields and using (ITT]) . (H9l) and (l20]h we see that 

Nu — f (| VT| 2 ) dz. (21) 

Jo 

On the other hand, by the definition of Nu and (l20]h testing (TT8l) with u and 
using (1171) and (fT9]h we obtain 

1 f 1 

Nu=—— / (\u\ 2 )dz+l. (22) 

Ra J o 

We are now able to estimate the Nusselt number. Letting i e (0,1) be 
an arbitrary number, we estimate 

1 r e 1 

Nu < - ( wT ) dz+ - 

f Jo ^ 

as in the previous section. Because \T\ < 1 by the maximum principle for 
the temperature, we obtain via Jensen’s inequality, 

r l \ 1/2 i 

Nu ~ \£ J 0 <U?2) d j + A 
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By the Nusselt number representation (|22j) . this yields 


Nu< 


(. RaNu) 1 / 2 

W 2 



This bound is optimized by £ ~ 1, so that 


Nu < Ra 


because Ra^> 1. 

In a final step, we like to relate the current bound on the Nusselt number 
to a bound on the viscous dissipation rate. In fact, we will show that £ < 
( Ra) 2 Nu , so that the above statement turns into 

£ < (Ra) 3 . 

Indeed, we start noting that Darcy’s law QT8J) together with the boundary 
conditions (fHJl) and (OH provides us with Neumann boundary conditions for 
the pressure: d z \ z=0 p = Ra and d z \ z=1 p = 0. In particular, differentiat¬ 
ing the horizontal velocity components with respect to z, d z v = —'V y d z p, 
multiplying by v, averaging and integrating by parts yields 


(v ■ d z v) = —{v ■ V y d z p) = ((V y • v)d z p). 

In particular, the above values for d z p and the horizontal periodicity imply 
that 

(K'4t>)L. al =0. (23) 

It thus follows via integration by parts 


(|Vu | 2 ) dz = (v ■ d z v)\[ = l — I (v-Av)dz^— I (v ■ Av) dz. (24) 


Now notice that —Av = (Ra)V y d z T. Indeed, taking the divergence of (fT8l) 
yields A p = ( Ra)d z T , and thus —Av = V y A p = (Ra)W y d z T. Therefore, 
becomes 


(| Vv| 2 ) dz = Ra / (v ■ V y d z T) dz — —Ra / ((V„ ■ v)d x T) dz. 


Since |V y • v\ < |Vyu| < |Vu|, we can use the Cauchy-Schwarz inequality to 
obtain 




1/2 


<| Vup) dz < Ra [ I (|Vn | 2 ) dz J (|VT| 2 ) dz j , 


and thus, via (1211) . 


(| Vu| 2 ) dz < (Ra) 2 Nu. 
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The estimate of the vertical velocity component is easier because of (1191) . 
Indeed, testing Aw = (Ra)AT — d z Ap = ( Ra)A y T with w, integrating by 
parts and using Young’s inequality yields 



Combining the last two estimates finally yields £ < ( Ra) 2 Nu as desired. 
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